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Abstract 

Parity (P) and time (T) invariance violating effects in the Ra atom are 
strongly enhanced due to close states of opposite parity, the large nuclear 
charge Z and the collective nature of P, T-odd nuclear moments. We have 
performed calculations of the atomic electric dipole moments (EDM) pro- 
duced by the electron EDM and the nuclear magnetic quadrupole and Schiff 
moments. We have also calculated the effects of parity non-conservation pro- 
duced by the nuclear anapole moment and the weak charge. Our results show 
that as a rule the values of these effects are much larger than those considered 
so far in other atoms (enhancement is up to 10 5 times). 
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I. INTRODUCTION 



The lower energy levels of radium corresponding to configurations of different parity have 
very close energies. This leads to a strong enhancement of the various parity (P) and time 
(T) invariance violating effects. In our previous paper Jl| we considered the states 7s6d 3 D 2 
with E = 13993. 97cm~ 1 and 7s7p 3 Pi with E = 13999. 38cm _1 , which are separated by 
a very small interval of ~5 cm -1 (~ 1CT 3 eV). Simple estimates showed that the effects 
of nuclear P- and T-odd moments such as the magnetic quadrupole moment (MQM), the 
Schiff moment (SM) and the anapole moment (AM) are many times larger than in all atomic 
systems considered before. In the present paper we present more accurate calculations of 
these and other parity and time invariance non-conserving effects in those states of the 
radium atom where the effects are large. We use a relatively simple ab initio approximation 
to perform the calculations. The approximation is a reasonable compromise between the 
simplicity of the calculations and the accuracy of the results. It is based on relativistic 
Hartree-Fock (RHF) and configuration interaction (CI) methods. A minimum number of 
basis states are used at the CI stage of the calculations. However, important many-body 
effects, such as polarization of the atomic core by an external field and correlations between 
core and valence electrons, are included in the calculations of single-electron matrix elements. 
To control the accuracy of the calculations we also calculated hyperfine structure intervals 
and lifetimes of lower states of radium and its lighter analog barium. 

Our calculations confirm the estimates done in the previous work [I| and show that the 
value of most P- and T-odd effects in radium is much higher than in other atoms considered 
before. The parity non-conserving (PNC) electric dipole transition amplitude between the 
ground and 3 Di even states is about E1 PNC ~ 0.8 x 10~ 9 (Q w /N)iea , which is 100 times 
larger than the measured PNC amplitude in cesium || and about 5 times larger than the 
corresponding amplitude in francium ||. The enhancement of the electron electric dipole 
moment (EDM) in the 3 Di state of Ra is about 5400, which is again many times larger than 
corresponding values for the ground states of Fr (910) and Au (260) 0. The transition 
amplitude between the ground and 3 D 2 even states induced by the nuclear anapole moment 
is about 10~ 9 eao, which is more than 10 3 times larger than a similar amplitude in Cs [Q. 
Also, the EDM of the Ra atom in the 3 D 2 state induced by the nuclear Schiff and magnetic 
quadrupole moments is strongly enhanced. Both contributions (SM and MQM) are about 
10 — 19 ?7 e • cm (rj is the dimensionless constant of the P-, T-odd nucleon-nucleon interaction). 
This is again about 10 5 times larger than the EDM of the Hg atom which currently gives the 
best limit on rj ||. All this makes radium a very promising candidate for the experimental 
study of P- and T-odd forces by means of atomic physics. 

II. METHOD 

We use relativistic Hartree-Fock (RHF) and configuration interaction (CI) methods to 
construct two-electron wave functions of the ground and lower excited states of barium and 
radium. The calculations start from the RHF method for a closed shell system corresponding 
to the ground state configuration (6s 2 for Ba and 7s 2 for Ra). Since nsnp and ns(n — l)d 
configurations, with n = 6 for Ba and n = 7 for Ra, do not correspond to a closed-shell 
system, we calculate p and d basis states in a model HF potential. For example, to calculate 



2 



7p and Qd states of Ra, we keep all other states frozen, remove the contribution of one 
7s-electron from the direct HF potential and use this potential to calculate the required 
states. The same procedure applies for Ba. Thus, we have five single-electron basis states 
for the CI calculations (nsi/2,npi/ 2 ,np 3 /2, (n — 1)^3/2, (n — l)d 5 / 2 ). It turns out, however, 
that this simple CI approximation significantly overestimates the relative value of spin-orbit 
intervals for the odd-parity states and underestimates it for the even-parity states. This 
affects the accuracy of the calculation of P- and T-odd effects because most of them involve 
transitions with a change of spin which are sensitive to the value of the relativistic effects. We 
found that the spin-orbit intervals are sensitive to the screening of the Coulomb interaction 
between two external electrons (recall that Coulomb integrals contribute to the spin-orbit 
splitting due to the difference between the single-particle radial wave functions belonging 
to different components of the single-particle doublets). To improve the quality of the wave 
functions we introduce fitting factors to the Coulomb interaction in the CI calculations (k 
is the multipolarity of the Coulomb interaction). It was found that multiplying all Coulomb 
integrals of multipolarities 0, 1 and 2 by factors f = 0.7, fx = 0.75, f 2 = 0.9 significantly 
improves the energies and fine structure intervals of lower odd and even states of barium 
and radium. These factors simulate the effect of the screening of the Coulomb interaction 
between valence electrons and core electrons. They also compensate to some extent the 
effect of the incompleteness of the basis set. 

To calculate values other than energy, such as the effect of electron interaction with 
photons and nuclear P- and T-odd fields, we also include core polarization effects (direct 
and exchange RPA-type corrections) and core-valence correlation effects (the Bruckner-type 
correlation corrections). These two effects are very important for the considered states of 
radium. Indeed, consider mixture of the 3 Dj and 3 Pj> states by the P- and T-odd interac- 
tion W. Corresponding dominant configurations (7s6d and 7s 7p) can only be mixed by a 
(7p\W\6d) matrix element. However, this matrix element is extremely small in the Hartree- 
Fock approximation. This is because the electron interaction with P- and T-odd nuclear 
moments is localized in the vicinity of the nucleus where the (i-electron does not penetrate 
due to the centrifugal barrier. On the other hand, the polarization of the electron core by 
these moments produces a long-range correction 5V to the HF potential which effectively 
renormalizes the interaction of an external electron with the nucleus. The corresponding 
matrix element {7p\W + 5V\6d) is not small even in the case of the p — d transition due 
to the long-range of the renormalized interaction W + SV. Note that (7s\W\7p) matrix 
elements also contribute to the mixture of the 3 Dj and 3 Pj> states due to the configuration 
interaction. Thus, there is an interference of several factors: the s — p matrix elements are 
large but their contribution is suppressed due to the smallness of the configuration mixing. 
The p — d matrix elements are considerably smaller (although not negligible) but they ap- 
pear in the dominating configurations. It cannot be said in advance which transitions are 
more important and as we see from our calculations there are cases when s — p transitions 
dominate over p — d and vice- versa (see below). The Bruckner- type correlation corrections 
(the correlation corrections to the single-electron wave functions) are also important, since 
they increase the density of an external electron on the nucleus by ~ 30% (see e.g. |J). 

The full scale inclusion of the core polarization and correlation effects into the CI cal- 
culations (see, e.g. 0) lies beyond the framework of this research. We adopted a simplified 
approach in which the corresponding corrections are calculated for the single-particle matrix 
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elements. The relative values of the renormalization of the matrix elements by the core po- 
larization and core-valence correlations have been extrapolated from accurate calculations 
of the core polarization and Bruckner-type correlation corrections for the radium positive 
ion. Ra + has a simple electronic structure - one electron above closed shells - and the 
corresponding procedures are well defined for it ||. 

To check our method and the accuracy of our results, we calculated the hyperfine struc- 
ture (hfs) constants of 213 Ra and 137 Ba. The results for the energies and hfs constants are 
presented in Table |[ One can see that even for this very simple CI approximation the 
accuracy of the energies and fine structure intervals is very good. The accuracy of the hfs 
constants is also good for the most important states 3 D 2 and 3 Pi. Table [IT] shows the effect 
of the core polarization (RPA) and Bruckner-type correlations (E) on the single-electron 
matrix elements. One can see that these effects play a crucial role in the p — d matrix 
elements. However, their contribution to the s — p matrix elements is also very important. 



III. PARITY VIOLATION IN 7S 2 -> 7S6D TRANSITIONS 



A. Spin-independent parity non-conservation 



The Hamiltonian Hpnc of the interaction of an electron with the nuclear weak charge Qw 
(formula ( |^6T ) in the appendix) mixes states of the same total momentum J and opposite 
parity. Thus, electric dipole transitions between states of initially equal parity become 
possible. In particular, the transition between the ground state 1 S and the excited 3 Dx 
state is enhanced due to the closeness of the opposite parity state 1 Pi. The dominating 
contribution to this transition is given by 

(7s 2 1 S \d z \7s7p 3 P 1 }(7s7p 3 P 1 \H PNC \7s6d 3 D 1 } 

ElpNC ~ Epnj-Eppj • (1) 

Apart from the enhancement, there are several suppression factors in (|I]). First, the electric 
dipole matrix element is small because of a change of spin. It is 3 to 5 times smaller than 
most of those amplitudes which do not change atomic spin (see Table fTJ). Second, in the 
matrix element of the PNC interaction, leading configurations produce only the p 3 / 2 — d 3 / 2 
single-electron matrix element which is small. It is not zero mostly due to core polarization. 
However, it is about 25 times smaller than the su2 — p\i% matrix element. The latter 
contribute to the PNC amplitude due to configuration mixing. Our calculations show that 
the contribution of the s — p transition to the PNC amplitude is about 7 times larger than 
the contribution of the p — d transition. In spite of some suppression, the final answer is 
quite large: 

225 Ra : El PNC = 0.77 x 1(T 9 (Q w /N)iea , (2) 
223 Ra : El PNC = 0.76 x 10~ 9 (Q w /N)iea . (3) 

This is one hundred times larger than the measured PNC amplitude in cesium and about 
5 times larger than the corresponding amplitude in francium ||. Even radium isotopes have 
close values of the amplitudes (approximately, the effect is proportional to the number of 
neutrons N). 
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B. Anapole moment 



The Hamiltonian of the electron interaction with the nuclear anapole moment is presented 
in the appendix (|A9|) . Similar to the spin-independent PNC interaction, it mixes states of 
opposite parity and leads to non-zero El-transition amplitudes between states of initially 
equal parity. However, it can also mix states with A J = 1 and it depends on the nuclear spin, 
so that its contribution to transitions between different hyperfme structure components are 
different. The corresponding expression is very similar to ([I]). However, dependence on the 
hyperfme structure must be included (see formula ( A12Q in the appendix for details). This 
amplitude is proportional to the ( 3 Pi | \ap(r) | | 3 D 2 ) matrix element. Contributions of different 
single-electron transitions into this matrix element are presented in Table [V[ Note the 
strong cancellation between terms corresponding to s — p and p — d transitions. This means 
that an accurate inclusion of the core polarization and core- valence correlation effects is very 
important indeed, as has been discussed above. We believe that the fitting of the energies 
helps to stabilize this matrix element similar to the case of the .El-transition amplitude. 

The results for 1 So— and 1 So— 3 D 2 transitions are presented in Table |V[ Note that 
the contribution of the anapole moment to the PNC amplitude ([3]) can be measured by 
comparing the amplitudes between different hyperfme structure components similar to what 
was done for cesium [fj| . However, it may be much more efficient to measure the effect of the 
anapole moment in the 1 So— 3 D 2 transition because it is about ten times larger due to the 
small energy denominator and because the nuclear spin independent PNC interaction does 
not contribute to this amplitude at all due to the large change of the total electron angular 
momentum A J = 2. 



IV. ATOMIC ELECTRIC DIPOLE MOMENTS 
A. Electron EDM 

An electron electric dipole moment interacting with an atomic field mixes states with 
the same total momentum J and opposite parity. As a result, an atomic EDM appears. The 
EDM of radium in the 3 Di state is strongly enhanced due to the closeness of the opposite 
parity state 3 Pi. In an approximation when only the mixture of the closest states is included, 
the EDM is given by 

n (7s6d 3 D!| - er\7s7p 3 P 1 )(7s7p ^Vj\Hedm[IsM 3 Dx) 

s(3 Dl ) - ^Pi) ■ [) 

Calculations using formulae from the appendix give the following result 

d = 53704- (5) 

Note that a very strong enhancement is caused by the small energy denominator E( 3 Di) — 
£( 3 Pi) = 0.001292 a.u. 
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B. Schiff moment 



Electron interaction with the nuclear Schiff moment also produces an atomic EDM. The 
EDM of Ra caused by Schiff moment is strongly enhanced in the 3 D 2 state. Its value is 
approximately given by 

, _ n (7s6d 3 B 2 \d z \7s7p 3 P 1 ){7s7p ^.{HsM^sGd 3 D 2 ) 

More detailed formula which include the dependence of (§) on the hyperfine structure is 
presented in the appendix QA19j) . 

Table |I^ shows single-electron contributions to the ( 3 P 1 ||if5j\,/|| 3 D 2 ) matrix element. 
Note that s — p transitions strongly dominate here. However, the contribution of the p — d 
transitions is not negligible and should be included for accurate results. 

Calculated values of the radium EDM induced by the Schiff moment are presented in 
Table IVTL 



C. Magnetic quadrupole moment 

Electron interaction with the nuclear MQM can also produce an EDM of an atom. 
However, in contrast with the case of the Schiff moment, the MQM of isotopes where the 
nuclear spin I < 1 (like 225 Ra, where I = 1/2) is zero. The EDM of Ra in the 3 D 2 state is 
given by a formula similar to (^) 

, _ n {7s6d 3 D 2 \d z \7s7p 3 P 1 )(7s7p 3 Pi|ff M qM|7s6rf 3 D 2 ) 
dz ~ Z £( 3 D 2 ) - i^) ■ [0 



Again, more detailed formula can be found in the appendix (|A24 ) . 



Table [TV] shows single-electron contributions to the ( Pi||A TO fc|| D 2 ) matrix element. Note 
that in contrast to the cases of the Schiff and anapole moments, p — d transitions dominate 
over s — p transitions in this matrix element. For the anapole moment these two types of 
transitions contribute almost equally, while for the Schiff moment s—p transitions dominate. 
Note that s — p transitions appear due to configuration mixing only, while contribution of 
the p — d transitions is extremely small if core polarization is not included. This indicates 
once more that even for a rough estimation of the time or parity invariance violating effects 
in Ra an inclusion of the appropriate many-body effects is essential. 

Calculated values of the radium EDM induced by the magnetic quadrupole moment are 



presented in Table VII 



V. LIFETIMES 

To plan experimental measurements of space and time invariance violation in radium it 
is important to know the lifetimes of the states of interest. Apart from that, comparison of 
the calculated and experimental lifetimes can serve as a good test of the method used for 
calculation of P- and T- invariance violation since the same dipole transition amplitudes 
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contribute in either case. As far as we know, none of the radium lifetimes have been measured 
so far. On the other hand, some experimental data is available for barium. Therefore, we 
calculated lifetimes of lower states of both atoms. Results for dipole transition amplitudes 



are presented in Table |T| and the corresponding lifetimes are in Table [VII] . 

For the purpose of the present work, the most important states of radium are 3 Pi, 3 Di and 
3 D 2 states. The decay rate of the 3 Pi state is strongly dominated by the 3 Pi - 1 So transition. 
Transitions to the 3 Di and 3 D2 states are suppressed due to small frequencies. The 3 Pi - 1 So 
dipole transition amplitude involves a change of the atomic spin and therefore is sensitive 
to the value of the relativistic effects. This makes the amplitude numerically unstable. 
This probably explains the poor agreement between different calculations (see Table |V11J| ). 
However, we believe that the fitting of the fine structure which we have done for Ba and 
Ra (see Section |II|) brings the amplitude close to the correct value. This is supported by 
similar calculations for barium. The 3 Pi - 1 S amplitude contributes 38% to the decay rate 
of the 3 P! state of barium. Good agreement between calculated and experimental lifetimes 
of this state (see Table |VIII|) means that all transition amplitudes, including the 3 Pi - 1 So 
amplitude, are calculated quite accurately. 

The lifetime of the 3 Di state of Ra is determined by the 3 Di - 3 Po transition which is 
numerically stable. The lifetime of this state calculated by us is in good agreement with the 
estimations done by Budker and DeMille 0. 

The 3 D 2 state of radium is a metastable state. It decays only via electric quadrupole 
(E2) transition to the ground state. Calculations similar to the electric dipole transitions 
show that the lifetime of this state in the absence of external fields is about 15 seconds. 
However, measurements of the atomic EDM involve placing the atoms in a strong electric 
field. It is important to know how the lifetimes of the 3 D 2 and 3 D! states of Ra are affected 
by this field. The electric field mixes states of different parity and AJ = 0, ±1. If only an 
admixture of the nearest state is taken into account, the amplitude which determines the 
decay rate of a 3 Dj state is given by 

= CSol^Pi^PiKpDj) , , 

E( 3 Pi) — E( 3 Dj) • U 

Where £ is the electric field. This leads to the following decay rates 

W( 3 D 2 ) = 0.21£ 2 (9) 
W( 3 Dx) = 0.25 x 10~ 4 S 2 (10) 

For an electric field of 10 kV/cm, the lifetime of the 3 D 2 state is 30 /is, while the lifetime 
of the 3 Di state is 240 ms. This latter result is in good agreement with estimations done 
by Budker and DeMille @. Note that the state 3 D 2 , with maximum or minimal possible 
projection of the total momentum on the direction of the electric field (M = ±2), cannot 
be mixed by this field with the 3 Pi state. Therefore its lifetime is much less affected. 



VI. CONCLUSION 

The radium atom turns out to be a very promising candidate for the study of time 
and space invariance violating effects. All such effects considered in this paper are strongly 
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enhanced due to the high value of the nuclear charge Z and the closeness of the opposite 
parity states of the atom. Moreover, the contribution of different mechanisms to the time and 
space invariance violating effects can be studied separately if measurements are performed 
for different states and different isotopes of the radium atom. For example, the atomic EDM 
induced by the electron EDM is strongly enhanced in the 3 Di state, while contributions of 
the nuclear Schiff and magnetic quadrupole moments are strongly enhanced in the 3 D2 state. 
On the other hand, the magnetic quadrupole moment is zero for isotopes with nuclear spin 
/ = 1/2, like 225 Ra, while the Schiff moment for these isotopes is not zero. 

Calculations of the space and time invariance violating effects in radium reveal the im- 
portance of relativistic and many-body effects. The accuracy achieved in the present work 
is probably 20-30 %. However, a further improvement in accuracy is possible if such a need 
arises from the progress in measurements. 

This work was supported by the Australian Research Council. 



APPENDIX A: WAVE FUNCTIONS AND MATRIX ELEMENTS 

1. Radium wave functions 

Two-electron wave functions of the ground ( 1 So) and three excited ( 3 P 1 , 3 D 1 and 3 D 2 ) 
states of radium used in this work for the calculation of space and time invariance violation 
have the following form 

|7s 2 J = 0,L = 0,M = 0) = 

-0.9757|7si _i7si i) - 0. 115017m _i7m i) - 

2 ' 2 2 ' 2 2 ' 2 2 ' 2 

-0.0752(|7p 3 _37ps s> - \7ps_i7p-a i»+ (Al) 

2' 2 2'2 2' 2 2'2 

+0.0658(16^3 _36cb a) - \6da _i6ds i» + 

2' 2 2 ' 2 2 ' 2 2 ' 2 

+0.0702(|6rfs _s6d5 s> - |6rfs _ 3 6ds s) + 16^5 _i6rfs i», 

2 ' 2 2 ' 2 2 ' 2 2 ' 2 2' 2 2 ' 2 

|7s7p J = 1,L = 1,M= 1) = 

-0.9010|7si iTpi i) - 0.3537|7si _i7p 3 s> + 0.2042|7si i7p 3 i) - 

2'2 2 >2 2 ' 2 2 ' 2 2'2 2 '2 

-0.0976|7pi _i6rf 3 a) + 0.0563|7pi i6d 3 i ) + 0.0512|7p 3 _i6d 3 3 ) - (A2) 

1 2> 2 2'2 2'2 2'2 2' 2 2'2 

— 0.0591|7p 3 i6ds i) + 0.0512|7p 3 s6ds _i) - 0.0018|7p 3 _ 3 6rfs 5) + 

2 ' 2 2 ' 2 2'2 2 ' 2 2 ' 2 2 ' 2 

+0.0014|7p 3 _i6rfs 3) - 0.0010|7p 3 i6ds 1) + 0.0006|7p 3 3 6d s 

2 ' 2 2 >2 2'2 2 ' 2 2 ' 2 2 ' 2 

\7s6d J = 1,L = 2,M= 1) = 

-0.8660|7si _iQdz z) + 0.5000|7si i6rf 3 i ) + 0.0002|6rf 3 _ 3 6ds s ) - (A3) 

2 > 2 2 ' 2 ' 1 2 ' 2 2 ' 2 ' 2 > 2 2 ' 2 v 

-0.0001|6rf 3 _i6rfs s> + 0.0001|6rf 3 i6rfs i) - 0.0001|6rf 3 3 6^5 _i) - 

2 ' 2 2 > 2 2 ' 2 2 ' 2 2'2 2 ' 2 

-0.0021|7pi _i7p 3 s) - 0.0012|7pi i7p 3 i), 

li 2' 2 2'2' 1 2'2 2'2' 

\7sM J = 2,L = 2,M = 2) = 

-0.8087|7si i6rf 3 s> - 0.5366|7si _i6ds B ) + 0.2400|7si i6rfs s) - 

2'2 2 ' 2 ' 2 ' 2 2 ' 2 ' 2'2 2 ' 2 ' 

-0.0084|6rf 3 i6rf 3 2 > - 0.005916^3 _i6ds b) + 0.005316^3 i6ds 3 ) - (A4) 

2 '2 2 '2 2 ' 2 2 ' 2 2'2 2 >2 

-0.003216^3 3 6ds i) - 0.006816^5 _i6ds b) + 0.0091|6rfs 16^5 3) + 

2 ' 2 2 ' 2 2 ' 2 2 ' 2 2 ' 2 2 ' 2 
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+0.0130|7j9i i7p 3 a) + 0.003817^3 i7p 3 3). 

2 > 2 2 > 2 2 ' 2 2 > 2 

We use the following form for the single-electron wave function 

~ r { iag(r)Q(v/r) jlm ) ' (A5) 
Here a = 1/137.036 is the fine structure constant, Cl(r/r)ji m = —{a ■ n)0(r/r)^ TO . 

2. Spin-independent weak interaction 

The Hamiltonian of the spin-independent weak interaction of an electron with the nucleus 
is given by |TU 



G 

Hp NC = -^=p(r)g H /7 5 , (A6) 

where G = 2.22255 x 10~ 14 a.u. is the Fermi constant, p is the nuclear density (/ pdV = 1), 
Q w « —N + Z(l — 4 sin 2 9yy) is the nuclear weak charge, and 75 is a Pauli matrix. The 
matrix element of (IA6I) with wave functions (IA5I) has a form 



(jllim 1 \Hp NC \j 2 l 2 m2) = -i^QwRpNC$j-ij2$hlJ>rnxmv ( A7 ) 

Rpnc = aj p{f\9i - 9\f?)dr - radial integral, 
l = 2j-l. 

However, it is often more convenient to express (|A8|) in a form 

Vilim^Hpxclhhm*) = (-1)^™ ( jj^ ° * J H)a^=Q w C PNC RpNc, (A8) 



Cpnc = y 2j'i + ]-fijii2$ij 2 — angular coefficient for the reduced matrix element. 



3. Anapole moment 

The Hamiltonian of the interaction of an electron with the nuclear anapole moment has 



the form 111 



H AM = ^^ r ^K Ka p(r), (A9) 

where / is the nuclear spin, K = (i+|)(— 1) /+ ^~', I is the orbital momentum of the outermost 
nucleon, K a is a dimensionless constant proportional to the strength of the nucleon-nucleon 
PNC interaction |T^]. The matrix elements of the Hamiltonian ( |A9| ) between the many- 
electron states of the atoms depend on the hyperfine structure (see, e.g. W[ 
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(IJ'F\H AM \IJF) = ^^(-1)^' {jj> F } 



x ^I(I + l)(2I+l)(J'\\ap(r)\\J), (A10) 
F = I + J, J — atomic total momentum. 



The electron part of the operator (jA9Q is ap{r). Its single-electron matrix elements over 
states (POl ) have a form 



0'iiimi|ap(r)|j 2 / 2 m2) = (-l)*-"* f Jl 1 32 ) (C 1AM R 1AM + C 2AM R 2AM ) } (All) 

V —mi q m 2 / 



1 

2 




k 


J2 


i 

2 


1 


1 
2 






J2 


i 

2 


1 



1«2' 



-RiAAf = -4:ira J g 1 f 2 dr, 
R2AM = -4tt« J f x g 2 dr. 

The dominating contribution to the ^-component of the parity non- conserving electric 
dipole transition amplitude between the 1 S and 3 D! states of Ra induced by the anapole 
moment is given by 

F 1 F' 



eipv = i-i) F ~ f [If o f (- 1 ) 4F ' +J+J ' +2/+1 v%^v / iSy x 









I 




\F 1 


'A 




J' 





x ,/(2F + 1)(2F> + 1) ! j ) 1 } > F ) ^»™ggg^ . (A12) 
Here F = I + J, / = min(F, F'). 

4. Electron EDM 

The Hamiltonian of the interaction of the electron EDM d e with the atomic electric 



field E has the form JlO] 



H ED M = -deP(E--E), (A13) 

where 



P=[0-l)- S= io«?J' E=-W(r). 

The atomic EDM induced by ( |A13| ) can be calculated as an average value of the operator 
of the dipole moment over states mixed by an operator similar to ( |A13|) 

H' EDM = -d e (P-l)(V-E). (A14) 

Its single-electron matrix elements have a form 
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(jihmi\H' EDM \j2hm2) = (-l)ii-^ ( ° ^ j d e C EDM R EDM , (A15) 



dV 
dr ■ 



Redm = 2a 2 J gi^g 2 dr. 



Note that the selection rules and the angular coefficients are the same as for the spin inde- 
pendent weak interaction (|A"6|), while the radial integrals are different. 

5. Schiff moment 

The Hamiltonian of the interaction of an electron with the nuclear Schiff moment has 



the form Hl3 



H SM = 4irS-Vp(r), (A16) 

S = S1/I,S is Schiff moment. Many-electron matrix elements of ( |A16| ) depend on the 
hyperfine structure similar to 



(IJ'F\H SM \IJF) = (-1)™' [j j> F } S^ I(I+1 f I+1) (J>\\4nVp(r)\\J). (A17) 

The electron part of the operator ( |A16|) is 47r Vp(r). Its single-electron matrix elements over 
states (|A5|) have the form 

{hkm 1 \^Vp{r)\j 2 l 2 m 2 ) = (-l)i™ ^ _^ J ^ J C SM R SM , (A18) 
C SM = (-1)* + V( 2 ii + !)( 2 ^ + 1) ( 1 | I) C(h + k + l), 



1, if x is even 
0, if x is odd 



£(*) = 

Rsm = -47r/(/!/ 2 + a 2 9l g 2 )f r dr. 



The EDM of Ra induced by the nuclear Schiff moment for a particular hyperfine structure 
component of the 3 D 2 state is approximately given by 

J' J F \ [ I I 1 



J _ 9 I - - | r_1 A2F+2/+J+J' J J 1 r [ ) 1 1 1 L l / (J+D(2J+l) 

(2F | l)5 (7a6d 3D j|l gl ll 7s7 P 3P ^)( 7g7 P 3 P.,/||47rVp(r)||7 S 6d 3 Dj) x^gv 



6. Magnetic quadrupole moment 

The Hamiltonian of the interaction of an electron with the nuclear magnetic quadrupole 



moment has the form 13 
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M 

AI{2~T- I) 



HmQM — — 7777T7 TrtjnkAmk, 



2 

tmk = Imh + him ~ -&k m I(I + 1), (A20) 

r 



Its many-electron matrix element is 

3M 



(IJ'F\H MQM \IJF) = ^L_^5(2F + 1)(2J + 3) (J + 1)(2J + 1)7(2/ - 1) x 

2 2 ] 

J' I F \(J'\\A mk \\J). (A21) 
J / F J 

(A22) 

The single-electron matrix elements of the operator A mk have the form 

(MmilAnkLfefe/^) = (-1)*-"* ( Jl 2 J2 J (Ci A/Q M + ^2A/QM)^AfQAf , (A23) 

\ —777,1 9 m 2 ) ' ' 

CiMQM = (-l)^-*| v /(2 Jl + 1X272 + 1) ( j | ^1+^2 + 1), 



Cimqm = (-l) J1+j2+ ' 2+1 4 v /5(2j 1 + l)(2j 2 + l)(2Zi + 1)(2Z 2 + 1) x 

[Ik k 



h 1 z 2 





2 ji J 2 
1 1 1 



2 2 

Rmqm =ocJ F(r)(g 1 f 2 + /i£ 2 )<ir, 

r/r^, if r < 
1/r 3 , if r > ' 

r^v — nuclear radius. 



where .F(r) 



The EDM of Ra induced by the nuclear MQM for a particular hyperfme structure component 
of the 3 Dj state is approximately given by 



4 = 2 ( F ^{;V-1)^(2F + 1)I^ 



(21 + 3) (7+1) (21 + 1) 

X 



— F F J K v ^ ; 4 ^ ^(2/-l) 

(7s6rf 3 Dj||£l||7 S 7p 3 Pj')(7s7p 3 Pj,p m fc||7s6d 3 Dj) 




For the case of the EDM in the 3 Dj state, J' = 1, J = 2 in flA2|). 



. (A24) 
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TABLES 



TABLE I. Energies and hyperfine structure constants of lower excited states of 137 Ba (1=3/2 
/Lt=0.937365) and 213 Ra (1=1/2,^=0.6133). 

Atom State Energies (cm -1 ) hfs constant A (MHz) 







Calc. 


Exper. 


Calc. 


Exper. 


6s5d 


3 D! 


9225 


9034 


-632 


-520.5 




3 D 2 


9346 


9216 


357 


415.9 




3 D 3 


9554 


9596 


504 


456.6 




X D 2 


12147 


11395 


-26 


-82.18 


6sQp 


3 Po 


12203 


12226 








3 Pi 


12577 


12637 


1233 


1150.59 




3 P 2 


12464 


13514 


878 






x Pi 


18042 


18060 


-48 


-109.2 


7 'sip 


3 Po 


12971 


13078 








3 P1 


13926 


13999 


8058 






3 P 2 


16660 


16689 


4637 






X P1 


21033 


20716 


-1648 


-2315 


7s6d 


3 Dx 


13893 


13716 


-4108 






3 D 2 


14042 


13994 


1749 






3 D 3 


14299 


14707 


2744 






X D 2 


17750 


17081 


-320 
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TABLE II. Single-electron matrix elements of the P- and T-odd interactions for radium (pre- 
sented reduced matrix elements of an electron part of the Hamiltonian as specified in the table, 
see Appendix for details). All values are in atomic units. 



Matrix element 




Approximation 


Even or 




RHF a 


RHF+RPA+S b 


odd c 



Spin-independent PNC interaction, H = p(r)^ 



(7si /ol 
\ 1/^1 


\H\ 


l 7 Pi /?) 

I " ± / Z / 


-2769 


-3832 


Odd 


( 7 P3/2 


\H\ 


6^3/2) 


0.004 


-146.8 


Odd 








Nuclear Anapole moment, H = ap(r) 




( 7s l/2 


\H\ 


7 Pl/2) 


-503 


-577 


Odd 


( 7s l/2 


\H\ 


7 P3/2) 


-0.508 


20.26 


Even 


( 7 Pl/2 


\H\ 


6^3/2) 


-0.024 


-66.29 


Even 


( 7 P3/9 

\ Jr Of Zr 1 


\H\ 


|6cL/o) 
1 0/ z / 


0.0006 


-29.99 


Odd 


( 7 P3/2 


\H\ 


6^5/2) 





11.71 


Even 








Electron dipole moment, H = 


(0- 1)SE 




( 7s l/2 


\H\ 


7 Pl/2) 


12.06 


17.05 


Even 




\H\ 


|6rf 3 /2> 


0.556 


2.082 


Even 








Nuclear Schiff moment, H = 


47rV / o(r) 




( 7s l/2l 


\H\ 


|7Pl/2> 


-44400 


-63027 


Even 


( 7s l/2l 


\H\ 


l 7 P3/2> 


-32550 


-56730 


Odd 


<7Pl/2 


\H\ 


1 64/2 > 


-1497 


1873 


Odd 


<7P3/2 


\H\ 


M3/2) 


-0.03 


2767 


Even 


<7p 3 /2 


\H\ 


|6d 5 / 2 > 


-0.07 


8163 


Even 








Nuclear Magnetic quadrupole moment, H = A m k 




( 7s l/2 


\H\ 


|7p 3 /2> 


17.28 


25.06 


Odd 




\H\ 


|6d 3 / 2 ) 


2.831 


2.933 


Odd 




\H\ 


1 6^5/2) 


-0.2017 


6.631 


Even 


(7^3/2 


\H\ 


|6d 5 / 2 > 


0.5389 


4.011 


Odd 



a Relativistic Hartree-Fock 

b Core polarization and core-valence correlation interaction are included 
c Even means that = (j||.£f||i); odd means that = — 
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TABLE III. El-transition amplitudes for Ba and Ra (|(2||4||i)| a o)- 



Transition 




Ba 




Ra 


i 


j 


A TTi"nlitnnp 


T^Ypm i pti pv f ■ — f ■ la 11 I 

±. L CVJ LI OX J- V Cj, Co I Gu* LL> / 


A TDTilitii ri p 


T^YpmiP'npv f- — f- (pi i 


3 Po 


3 Di 


2.3121 


n m 473 


3 D449 


-D Df)29D4 


3 Pi 


% 


D 4-537 


D 05758 


1 D337 


D D6379 


3 Pi 


3 Di 


2 D1 D8 


n m 64i 


2 6389 


n nm 299 


3 Pi 


3 D 2 


3.4425 


n m 55Q 


4 4399 


n nnnn947 


3 Pi 


X D 2 


o i fii n 


D DD566 


D D467 




3 P 2 


3 Di 


0.5275 


02042 


0.7166 


01354 


3 P 2 


3 D 2 


9 D94 






01 998 


3 P 2 


3 D 3 


4.777 


0.01785 


6.3728 


0.009027 


3 P 2 


X D 2 


0.1573 


0.00966 


0.1499 


-0.001790 


x Pi 


x So 


5.236 


0.08229 


5.4797 


0.09439 


x Pi 


3 Di 


0.1047 


0.04113 


0.4441 


0.03189 


x Pi 


3 D 2 


0.4827 


0.04030 


1.188 


0.03063 


x Pi 


X D 2 


1.047 


0.03037 


2.4053 


0.01656 



TABLE IV. Single-electron contributions to the two-electron matrix element ( 3 Pi||JY|| 3 D 2 ). 
Dash means no contribution due to selection rules. Zero means very small contribution. Same 
units as in Table [n| 



Transition 


H = -er a 


H = ap(r) h 


H = 4vrV/3(r) c 


H — A m ]^ 


7s i/ 2 - 7pi/ 2 


-0.3677 


58.78 


6421 




7pi /2 - 7s i/ 2 


0.0215 


3.431 


-375 




7s 1/2 - 7p 3/2 


-0.1306 


-0.515 


1441 


-1.565 


7p 3 / 2 - 7s 1/2 


-0.0585 


0.230 


645 


-0.289 


lp l/2 - 6d 3/2 


0.0020 


0.029 


-1 


-0.003 


64/2 - 7pi/ 2 


3.856 


-54.08 


-1528 


-1.848 


7pi/ 2 - 64/2 











64/2 - 7p 1/2 








-2.032 


7p 3/2 - 64/2 


0.0004 


0.006 


-1 




64/2 - 7p 3 / 2 


-0.2470 


3.522 


325 




7p 3/2 - 64/2 














64/2 - 7p 3 /2 


1.364 


2.442 


-1702 


-0.736 


Total 


4.4399 


13.85 


5226 


-6.473 



a For El transition amplitude 

b For anapole moment contribution 

c For Schiff moment contribution 

d For Magnetic quadrupole moment contribution 
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TABLE V. Parity non-conserving El-transition amplitude induced by nuclear anapole moment 



I 


F 


F' 


(dz) 


in units 10 10 K a ieao 












0.5 


0.5 


1.5 


2.05 


-20.3 


1.5 


1.5 


0.5 


-0.58 


5.7 




1.5 


1.5 


-1.4 


13.8 




1.5 


2.5 


1.3 


-12.9 



TABLE VI. EDM of Ra atom in the 3 D2 state induced by nuclear Schiff moment 





F 




d z (a.u.) 




d z (ecm) 




0.5 


1.5 


-0.94 x W 8 S 


-0.19 x 10" 


-ll^a 


-0.36 x 10" 




1.5 


0.5 


-0.16 x 10 s S 


-0.42 x 10" 


-ll„b 


-0.80 x 10- 


-19 

7] 


1.5 


1.5 


-0.30 x lO 9 ^ 


-0.81 x 10" 


-ll^b 


-0.15 x 10- 


- 18 7] 


1.5 


2.5 


-0.28 x 10 9 S 


-0.76 x 10- 


-ll^b 


-0.14 x 10- 


- 18 7] 



a Nuclear Schiff moment S is assumed to be S = 400 x 10 8 ?7 e fm [ 20 1 
h S = 300 x 10 8 ?? e fm 3 @ 



TABLE VII. EDM of 223 Ra isotope (7 = 3/2) in the 3 D 2 state induced by nuclear magnetic 
quadrupole moment 



F 


4 a 


d z h 


0.5 


1344Mm e 


7.4 x W~ 20 r] e- cm 


1.5 


1292Mm e 


7.0 x 10 _20 r? e- cm 


2.5 


-806Mm e 


—4.4 x 10 _20 r/ e- cm 



a In terms of nuclear magnetic quadrupole moment M 

b M is assumed to be M = 10 _19 (r//m p ) e- cm, [21] where m p is the proton mass 
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TABLE VIII. Lifetimes of lower short-living states of Ba and Ra 



Atom 


State 


Lower states to decay to 




Lifetime 






via El-transitions 


This work 


Other data 


Ba 


3 Po 


3 Dr 


2.83 /us 






3 Pi 




1.37 /us 


1.2 /us a 




3 P 2 


SDi^Da^Dg/Da 


1.41 fj,s 






Pi 


lc 3^ 3n in 
L>1, L>2, L>2 


y.l ns 


8.5 ns a 


Ra 


3 Pl 


^o^Dx^Da 


505 ns 


420 ns b , 250 ns c 




3 P 2 


3 D!, 3 D 2 , 3 D3 


74.6 ns 






3 D! 


3p 

Po 


617 /us 


800 /us d 






3 Pl, 3 P 2 


38 ms 






X Pi 


iSo^Dx^Da^Da 


5.5 ns 





a Reference [1£] 
b Reference [17 



^Reference [18 



d Estimation, Reference 
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